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ABSTRACT 

Measurements of the ultrasonic attenuation of 10 and 30 Mcs/sec 
longitudinal waves by pulsed-echo techniques were made on a 99.999% 
pure single crystal of superconducting zinc in the [ 0001 | direction 
as a function of temperature from 4.2°K to 0.320°K using an open-ended 
type Hee cryostat. The attenuation was found to be frequency dependent 
and decreased less sharply near the superconducting transition tem- 
perature, T.> than predicted by the Bardeen-Cooper-Schrieffer (BCS) 
theory. Attenuation due to electron-phonon interactions only was 
found by subtracting from experimental points the value of residual 
attenuation gotten by extrapolation of the data to T = O°K. Using 
the BCS theory the Zero degree superconducting energy gap was found 


to be 2 €(0) = (3.36 + 0.13)kT, with T, = 0.817°K. 
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INTRODUCTION 


The rapid decrease in the attenuation of ultrasonic waves in 
superconductors below the superconducting transition temperature, We ; 


‘ Bommel investigated 


was first observed by Bommel and Mackinnon. 
the ultrasonic attenuation of shear and longitudinal waves in a cylin- 
drical crystal of lead in the frequency range from 9-27 Mcs/sec from 
approximately 18°K to 1.5°K. He found a marked decrease in the atten- 
uation in the superconducting state which was frequency dependent. No 
theoretical explanation was made, though he suggested that the results 
might lend support to the idea that superconductivity is due mainly to 
interactions between electrons and lattice vibrations. Mackinnon in- 
vestigated the absorption of 10 Mcs/sec sound pulses in a superconducting 
polycrystalline tin rod and in a polycrystalline magnesium rod between 
1.5°K and 4.2°K. In tin he found the attenuation decreased by a factor 


of about 50 and in magnesium that the attenuation decreased steadily 


with decreasing temperature. 


This decrease in the electron contribution to the ultrasonic atten- 
uation was first satisfactorily explained in 1957 by the theory of 
Bardeen, Cooper, and Schrieffer (BCS Leer. The BCS theory predicts 
that the ratio of the electron contribution to the attenuation in the 


superconducting state, q. , to that in the normal state, q mas 


- Zz 


Q e(T : 
n exp liar | =e il 
where 2€(T) is the temperature dependent energy gap. The theory further 


predicts that as T approaches zero, the energy gap approaches 3.50 kT. 


4 : : 
Morse and Bohm were the first to apply this theory to the analysis 
of ultrasonic attenuation data on tin and indium, shortly before the 


publication of the BCS theory, basing their analysis on a conversation 
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with Professor Bardeen. Their analysis pointed out the remarkable 
agreement between the BCS theory and their experimental data. Since 
then the BCS theory has been applied to many superconducing Healer ae 
These investigations have all shown close agreement between the 
theory and experiment, but have also indicated the need for refinement 
of some of the simplifying assumptions of the BCS theory. The theory 
shows that the basic interaction causing superconductivity is that of 
the exchange of a virtual phonon between two electrons. An electron 
moving through the lattice creates a virtual phonon which is absorbed 
by a nearby electron. If the distortion of the lattice is such as to 
surround the electron by a positive charge greater than the electronic 
one, the second electron sees an attractive potential greater than the 
repulsive Coulomb one and is attracted. These electrons form what is 
known as Cooper pairs. The theory then shows that the lowest energy of 
the system is achieved by the maximum pairing of electrons of opposite 
spin such that their momenta add to zero. The difference in the normal 
ground state energy and the Superconducting ground state energy is due 
entirely to an interaction potential, less than zero, arising from the 
correlation between the Cooper pairs. The energy of this interaction 
leading to the scattering of one pair to another is the difference 
between the energy of interaction by means of a phonon and the energy 


of Coulomb interaction. 


The existence of an energy gap then follows from the small, but 
finite, amount of energy required to break up a Cooper pair. But more 
important, having a single unpaired electron removes a large number 
of pairs which could have occupied the state of a Cooper pair. This 
gap will be temperature dependent since increasing numbers of pairs 
will be broken up by thermal excitation as the temperature rises from 


O°K to 1. 
c 


Morse, Olsen, and Ga yancel pointed out the need to consider the 


anisotropy of the energy gap and showed that in tin there is considerable 


oe 





anisotropy (3.2 kT. a: Coe kT). In order to conserve both momentum and 
energy in a phonon-electron interaction, the component of the electronic 
velocity parallel to the direction of sound must equal the phonon veloc- 
ity; i.e., the sound velocity. Since the velocity of electrons at the 
Fermi surface is in general much greater than that of sound, only those 
electrons moving at essentially right angles to the direction of sound 
propagation can interact with the phonons. On the Fermi surface these 
electrons lie in a narrow band perpendicular to the sound velocity. 
Since the ultrasonic attenuation in a particular direction is an average 
over this narrow band, the energy gap measured will be an average for 
the electrons lying within the band. Anisotropies in the Fermi surface 
will then lead to anisotropies in the energy gap. The BCS theory makes 
the simplifying assumption that the interaction energy is isotropic and 
constant near the Fermi surface. It is this assumption that leads to 
the BCS prediction that the zero degree energy gap is the same for all 


superconductors and is equal to 3.5 kT: 


The relation for the ratio of ultrasonic attenuation in the Super- 
conducting state to that in the normal state follows from the coherence 
of paired electrons. In the normal state scattering of an electron to 
another state is independent of any other transitions. However, the 
pairing of electrons in the superconducting state makes these transitions 
interdependent and the theory shows that the scattering by pairs is 
destructive. Thus, the absorption of phonons, causing ultrasonic atten- 
uation, decreases sharply below To: With the assumption of low energy 


eienons; i.e., Aw << 26€(0) , the ratio ae hoes follows. 


The present work on Zinc, which has a hexagonal crystal structure, 
was undertaken to determine the energy gap in the direction of the 
c-axis and to compare the results with measurements made in the base 
plane. Ultrasonic frequencies on the order of 10/ cps were used, 


well satisfying the condition for the given expression of of fo A 





EXPERIMENTAL PROCEDURES 


A cylindrical, single crystal of zinc, approximately 2%" long and 
3/4" in diameter with the c-axis, [0001] direction, parallel to the 
cylinder axis was purchased from Semi-Elements, Inc. The purity of the 
sample was stated to be in excess of 99.999%. From this ingot the experi- 
mental sample was prepared to the approximate desired size, limited by 
the cryostat, by cleaving both ends and etching. The ingot was held in 
a styrofoam form, immersed in liquid nitrogen and cleaved by laying a 
new, sharp, razor blade at the desired point and rapping it with a 
mallet. To reduce the diameter, the ends were coated with brushing 


lacquer and the ingot etched in dilute HCl acid. 


Standard ultrasonic pulse-echo Pacihmdantes: were used and will be 
discussed only in terms of the experiment. The success of this method 
depends strongly on having the end faces of the sample flat, parallel, 
and perpendicular to the direction in the crystal one wishes to investi- 
gate. To achieve this, the zinc sample was cast with a cerro-alloy in 
a bored out, cylindrical jig whose faces were machined flat, parallel 
and perpendicular to the crystal axis to a tolerance of 0.003". The 
ends were first lapped on metallographic paper and then polished on a 


diamond dust wheel. 


The final sample was (3.1602 + 0.0002) cm long and approximately 
5/8'' in diameter at room temperature. A Laue back reflection pattern 
of the zinc crystal showed concentric circles, indicating a thin poly- 
crystalline layer on each end. However, the x-ray spots showing through 


the circles indicated that the faces were perpendicular to the c-axis. 


A 4", X-cut, 10 megacycle quartz transducer was next bonded to the 
zinc sample. The success of this bond is a critical step in the experi- 
mental procedure. After some experimenting, the following procedure 
proved to be reliable. Both the quartz and sample faces are cleaned first 


with acetone to remove greasy substances, second with alcohol to remove 


mre 





water, and finally with acetone to ensure the surfaces are dry. The 
best bonds were obtained on days when the humidity was low. Next a 

very small amount of analytical grade glycerin was applied to the zinc 
and the transducer centered and squeezed to make the glycerin flow and 
also to remove air bubbles. Glycerin proved to be the most satisfactory 
bonding material, though its use required the making of a new bond for 
each run. The critical temperature for the bond was from 40°K to 50°K. 
If the bond held on being cooled through this range; the bond was good 
for the run. However, as can be seen, one has already started to cool 
with oe and if the run must be stopped for a broken bond, the delay is 


both costly and time consuming. 


A schematic for the ultrasonic pulsing equipment is shown in Fig. 1. 
Radio frequency pulses (10 Mcs/sec and 30 Mcs/sec) of approximately two 
microseconds' duration were generated by an Arenberg Pulsed Oscillator 
triggered by a Time Mark Generator with a repetition rate of 1000 cps 
or 100 cps. The slower rate had to be used below 0.5°K to reduce heating 
of the sample. The rf electrical pulses were simultaneously fed to a 
Type RM 45A Tektronix Oscilloscope, externally triggered by the Time 
Mark Generator, and impressed across the piezoelectric quartz transducer. 
The zinc furnished the rf ground and the oscillator, tuned to the funda- 
mental or odd harmonic of the transducer, caused it to vibrate and send 
sound waves through the zinc. The sound waves reflected off the other 
end return to the transducer to produce a voltage across the electrodes 
furnishing the initial pulse and also to be reflected back through the 
sample. With a pulse rate of 1000 cps, these reflected echoes had 
adequate time to decay away between pulses. The voltage produced by each 
returning echo was amplified and displayed in proper time relation to 
the initial pulse on the oscilloscope. The change in height of each echo 
is a measure of the attenuation in the sample due to phonon-electron 
interaction, impurities, imperfections and in the transducer bonding 
system. Absolute attenuations are gotten by taking pictures of the scope 


traces with a Polaroid Oscilloscope Camera, Model 196A, Special, using 
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3000 speed/Type 47A Polaroid film and plotting echo heights vs echo 
number on semi-log paper. Attenuation effects other than those due to 
phonon-electron interactions are subtracted off as discussed in the 


treatment of the data. 


Data were taken in the temperature range from 4.2°K to approximately 
0.320°K. The temperature range was obtained in a ee cryostat of the 
open-ended type, a schematic of which is shown in Fig. 2. It is desirable 
to cool the sample as fast as possible to reduce cooling time, but this 


is tempered by the necessity of not straining the bond. 


The sample was loaded in a copper sample holder, inserted in the 
Hee cryostat and the Hee chamber flushed with dry nitrogen and left at 
one atmosphere of dry nitrogen. The vacuum jackets of the lee cryostat 
and the He" dewar were thoroughly flushed with dry nitrogen and left 
soft at approximately 300 microns dry nitrogen pressure. The si dewar 
was also flushed and left at one atmosphere dry nitrogen pressure. The 
nitrogen dewar was then filled with liquid nitrogen and the whole system 


cooled to approximately 77°K which takes about six hours. 


The Hee and ne jackets were then pumped out until the pressure was 
less than 5 x NO om Hg. Keeping the pressure less than this in the Hee 
jacket was essential to prevent thermal conduction at the lowest temper- 
atures. While this was being done, the nitrogen and ee traps were 
filled and the Hed gas circulated for approximately one-half hour to 
condense out all impurities except Hees This was found to be necessary 
as impurities would freeze out in the cryostat and around the sample 
holder causing throttling when pumping full out on the ee The Hee 


impurity caused no trouble as thermometry was done using a carbon resistor 


3 
and not by using the He vapor pressure. 


When the jackets were pumped out and the Hee purified, the Hee dewar 
was filled and about eight mm Hg of fee let into the cryostat; the re- 
mainder being returned to the reservoir. The system was then cooled to 


M2 °K, taking about 12 hours. The process could be accelerated by 
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: 3 ; 4 eee ; 
softening the He jacket with He . However, achieving a satisfactory 
vacuum after the system was cooled down was very difficult since helium 
at that temperature is very difficult to pump. Without a good vacuum 


the low temperatures desired could not be reached. 


Once 4.2°K was reached, the Hee was pumped on and the Hee let 
into the cryostat. A minimum temperature of 1.09°K can be reached 
this way, well below the temperature at which Hee at one atmosphere 
will condense; i.e., 3.195°K. When equilibrium was reached, the cir- 
culating pump was valved off to pump on the Hees At 500 micons’ the 


booster pump was cut in and pumping full out on the fee the lowest 


temperatures were reached. 


Temperatures were determined by measuring the resistance of a 
%-watt Speer Carbon Resistor, nominal resistance 470 ohms, 10%, grade 
1002, in a Wheatstone Bridge Circuit employing a lock-in amplifier to 
measure the bridge unbalance. The resistor was calibrated against the 
vapor pressure of Hee down to 1.1°K using the "1958 We Scale of Temper- 
ature''. Down to the lowest temperature reached, the resistor was cali- 
brated against the magnetic susceptibility of cerium magnesium double 
nitrate (CMN), Ce,Mg,(NO,)).°24H,0 using a Cryotonics Mutual Inductance 
Pridge.° The susceptibility of CMN is strongly anisotropic with the 
susceptibility parallel to the symmetry axis essentially zero. To obviate 
the necessity of packing single crystals, powdered CMN was used. The 
effective susceptibility is then approximately 65% of the value perpen- 
dicular to the symmetry axis, but the gain in the packing fraction off- 


sets this decrease. 


Since the temperature dependent part of the susceptibility of CMN 


perpendicular to the symmetry axis follows Curie's Law 


X 2 


the mutual inductance bridge readings are correlated to the temperature 


Rilo 


and thereby calibrating the carbon resistor. 
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EXPERIMENTAL RESULTS 


Measurements of the ultrasonic attenuation of 10 and 30 Mcs/sec 
longitudinal waves were made in the [0001] direction on a single crystal 
of zinc whose ends were slightly polycrystalline. Data were taken in 
the temperature range heats approximately 0.320°K and in the earth's 
magnetic field. The superconducting transition temperature was determined 
to be T. = 0.817°K and the zero degree superconducting energy gap was 


measured as (3.36 + 0.13)kT.. 


To meaSure the attenuation as a function of temperature, pictures 
of echo trains were taken at different temperatures and the heights of 
the echoes plotted on semi-log paper to obtain the attenuation in units 
of attenuation per echo. These were plotted as a function of temperature 
and when extrapolated to T = o°K gave an estimate of the residual atten- 
uation. Subtracting the residual attenuation from experimental points 
gave the attenuation due solely to electron-phonon interactions. These 
data were then analyzed using the BCS function for reduced attenuation 


vs reduced temperature 


2 
exp ka el 


Vk = 


kT 


The details of all these analyses are given in Appendices I - III. 
In order to illustrate the methods, Figs. 3-14 graphically display the 


results for one set of data; those of 13 November 1964. 


After the transition temperature was determined and an estimate of 
the residual attenuation made, computer program RESATT was used to 
determine the residual attenuation and the zero degree energy gap. 
Figs. 3-5 show the point plots of reduced attenuation vs reduced 
temperature for three slightly different values of residual attenuation 
subtracted from the experimental data. By inverting the BCS expression 


the following function Y is defined: 
Y = 1/in(2/a* - 1) = (kT /e(t))t. 
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Figs. 6-8 depict the point plots of the function Y vs the reduced 
temperature, t , corresponding to the same values of residual atten- 
uation. All these figures have the corresponding curve plots for the 


BCS expression for 2¢€(0) = 3.50 kT. 


The criterion used to find the residual attenuation and the zero 
degree gap was that on the plot of Y vs t the curve best fitting the 
experimental points should become a straight line passing through the 
origin for t less then 0.35. This follows from the expression for Y 
which becomes linear for 2/a@* much greater than 1.00. The results of 


these analyses are given in Table I. 


Once the residual attenuation was determined, computer program 
VARIGAP was used to show how changing the zero degree energy gap would 
change the BCS curves for the residual attenuation and Y. Figs. 9-11 
show the BCS curve for reduced attenuation vs reduced temperature for 
a given residual attenuation and with changes of 0.10 kT. in the energy 
gap determined from RESATT. Figs. 12-14 give the corresponding results 
pore Y vs t. As can be seen the BCS curves do not shift much for the 


value of incremental change used. 


It was found that the attenuation was constant above Lye and frequency 
dependent. Table I shows that attenuation increased with increasing 
frequency. The higher attenuation for 30 Mcs/sec in the normal state 
on 5 January 1965 compared to that of 13 November 1964 reflects the 
greater attenuation losses from the poorer bond. The value for the 
zero degree gap was gotten from computations based on semi-log plots. 

It is felt that this method is better than that of comparing echo 
heights since the latter method involves propagating through all the 
data any error made in measuring the echo height chosen as the basis for 
comparison. However, Table I shows that essentially the same results 
are obtained using either method. The difference in the methods is 


apparent in the different values of residual attenuation obtained. Taking 
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the value at the lowest temperature reached for comparing echo heights 
required that the residual attenuation be determined from this value 
rather than from an average of low temperature values as was done 


using semi-log plots. 
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Frequency Normal Residual 2€(0) 
Mcs/sec Attenuation Attenuation KT 2€(0) 
a db/em a ,db/cm c ,,db/em =e 
n O é 
0.310 3.46 
30 3.092 0.320 3.58 
0,330 37,62 
0.525 2.96 
30 32350 0.565 3226 
0.605 3.49 
0.660 Pea fe 
10 1.376 0.700 3a25 
0.740 3.63 


Run 


Jan 


Jan 


Run 


Jan 


Jan 


65 


65 


65 


65 


TABLE 1 


RESULTS OF RESATT COMPUTATIONS 


Attenuation from Semi-log Plots 


Attenuation from Comparison of Echo Heights 


Frequency Normal Residual 2€(0) 
Mcs/sec Attenuation Attenuation KT 
Qa db/cm a ,db/cm ss 
n O 
0.290 3745 
30 3.092 0.300 3.58 
0.310 3265 
Oe 525 PEG es 
30 3.350 02505 2.99 
0.605 Sun oy 
0.700 Sa. 
10 1.376 Ones 55 
03730 


S058 


Results 


0.320 


0.565 


0.700 


Results 
a ,db/cm 
O 


0.300 


0.585 


Oey i> 


2€(0) 

KT 
c 

Cys: 


3.30 


ress) 





DISCUSSION OF RESULTS 


The value obtained for the transition temperature, Le = 0.817°K, 
is lower than those of other investigators. However, those other 
experiments were ceterminations of the zero magnetic field transition 
temperature for zinc and show considerable variation. In heat capacity 
measurements, PRilitns deeaeraned T. = 0.825°K. Geballe and Matthias!° 
in work on isotope effect measured T. = 0.855°K and Cochran and 
feeerher — in critical field experiments found T. = 0.875°K. If a 


parabolic relationship between reduced magnetic field and reduced 


temperature is assumed 


a ee a : 
H T : 
Cc Cc 


where He is the critical magnetic field above which zinc will not 
become superconducting, then the transition temperature is reduced 
about five millidegrees for He = 50 gauss and H = 0.5 gauss. Thus 


I = 0.817°K agrees favorably with the lowest value quoted above. 


Table II gives the results of several investigators who have 
measured the zero degree energy gap in Zinc. There have been no ultra- 
sonic attenuation measurements along the c-axis other than this one. 
The result of this work is higher than the value found by Biondi et al 
in a microwave absorption experiment in which there were ''violent 
oscillations" in the data for hv /kT less than 2.7. The measurement 
by Zemon and Boorse can only be taken as an average over all directions 
in the crystal. Bohm and Horowitz found the attenuation so great 
along the c-axis that they were unable to make any measurements. How- 
ever, their measurement of the energy gap in two directions in the base 
plane and the measurements of the energy gap along the c-axis by 
Biondi and this work show anisotropies in the energy gap. These results 
also indicate that the gap along the c-axis is less than that predicted 
by the BCS theory and also is less than that for directions in the 


base plane. 


Ge 








ZERO 


2 (0)/kT. 


3.8 + 0.2 
3.4 + 0.2 
3.00 + 0.15 
beet + 0.09 


peeoen 0.13 


TABLE Il 
DEGREE ENERGY GAP IN SUPERCONDUCTING ZINC 


Direction 


| 1270] 

| 1070] 

[0001 | 
Polycrystalline 


Sample 


| 0001) 





Ooo 


0.85 


0.838 


02625 


0.817 


Pye pe 


Method Investigator 
Ultrasonic Bohm and,, 
Attenuation Horowitz 


233 Mcs/sec 


Ultrasonic Bohm and 
Attenuation Horowitz 
233 Mcs/sec 


Microwave Biond}., 
Absorption et al 
15-100 Ges/sec 


Microwave Zemon d 
Absorption Boorse 
50-75 Gcs/sec 


Ultrasonic Present 
Attenuation Work 
10 & 30 Mcs/sec 





A more precise determination could be made if lower temperatures 
could be reached. The lowest temperature reached in this work was 
0.320°K; i.e., 0.38 t. Thus to determine the energy gap, the experi- 
mental data had to be extrapolated some 40% on the reduced temperature 


scale. 


The graphs of reduced attenuation vs reduced temperature can be 
fitted with a BCS curve quite well for low t. However, for t approaching 
1.00, the experimental points lie above the BCS curve. Bohm and Horowitz 
don't report the comparison of their points with the BCS curve. How- 
ever, investigators of other elements have found just the opposite; 
i.e., points for higher t lie below the BCS curve. For Niobiun, 
reper” found the effect to be frequency dependent. With Ty = 9.15°K, 
the BCS curve for the measured value of the energy gap, 3.63 kT. ; | 
fit the experimental points for 220 Mcs/sec longitudinal wave along 
the [100] axis. However, for 30 Mcs/sec the points lay below the BCS 
curve. He also found an isotropic result for the energy gap in the 
three crystallographic directions and suggested that any anisotropy in 
the gap may have been smeared out by impurity scattering. For thalliun, 
a = 2 38°K, Saunders and Paeonies found that attenuation of longitudinal 
waves decreased more sharply below T. than predicted by the BCS theory 


and that the anisotropy in the energy gap ranged from (3.76 - 4.10)kT . 


These and similar results on other elements are in contrast to the 
findings of this work which should be investigated further. One basic 
assumption tacit in all these computations was that the residual atten- 
uation was temperature independent. Though measurement above the 
transition temperature bore this out, the attenuation in the normal 
state for temperatures less than the zero magnetic field T. should be 
measured. This might be done by placing the sample in a magnetic 


field to destroy the superconducting state. 


Finally, measurements at higher frequencies should be carried out. 


. 1 ; ; at : 
Bobetic j has carried out the calculations giving the ration of @ /e 
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according to the BCS theory for an isotropic model at arbitrary 
frequencies and temperatures and predicted a discontinuity in attenu- 
ation at frequencies corresponding to the energy gap. This would 


: Z 
require frequencies of the order 10 Gcs/sec. 


In this work attempts were made to make measurements at 50 and 
70 Mcs/sec. However, so much power was needed to get an echo, that 
the sample heated excessively. Using a shorter sample would help 


solve this problem. 
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APPENDIX I 
DATA ANALYSIS 


Analysis begins by measuring the echo heights on each picture. 
This is most accuratelydone on an X-ray film reader to + 0.01 mm. These 
heights are then plotted on semi-log graph paper vs echo number; one 
graph per picture. The slope of this graph is absolute attenuation per 
echo. 

The basic problem in determining the absolute attenuation in the 
temperature range attained (4,2°K to approximately 0.320°K) is that the 
attenuation is so great from slightly below the superconducting transi- 
tion temperature, To to 4.2°K that one usually has only one echo. 
Below To» the attenuation drops rapidly and enough echoes appear to 
obtain the attenuation graphically. 

A typical plot of absolute attenuation, Q@, vs temperature, T, is 


as follows: 


| 

Ack : 

4 | 

Sete aieei oS = | 
a 

| 

X . 

| 
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The 
db/cm is 


We a 


Let 


absolute attenuation 


- residual attenuation 


attenuation due to conduction electrons in the superconducting 
state 


attenuation due to conduction electrons in the normal state 


: attenuation determined from raw data 


O 
Temperature - K 
Superconducting transition temperature 


Change in attenuation between the superconducting and normal states 


slope of the semi-log plots are in ee The attenuation in 
given as follows: 


Attenuation (db/em) = Slope of Semi-log Plot x 8.6858 


~ 2 x sample length in cm 


ssume an exponential decay of echo heights. 


A= Me 


amplitude of nha echo in the superconducting state 
amplitude of oi echo in the normal state 
amplitude of the initial pulse 

echo number 


: ; oO 
attenuation in the normal state from T. to 4.2 K: a constant 
as verified experimentally 


temperature dependent attenuation in the superconducting state 


sample length in cm 


-OT “n 
A‘e S 





a 
or om (db/cm) = = eee Ln (x aos + Ot; 


ae N 





To determine O,7 from the single echo pictures above TO: one com- 
Cc 
putes O,7 » with n= 1, for as many point by the above equation as one 


Cc 
can from the semi-log plots. The best horizontal straight line fit to 


a plot of these points gives Q@_. Then using the equation 


Te 


8.6858 Ay 


a , (db/cm ) = Aaa n ( Kae + 


N 


gives the attenuation for the remaining points below T. where there 
are not enough echoes to determine the attenuation from the semi-log 
plots. 

Extrapolating the plot of O. and 7 vst to T ="OMwriieriee an 


Cc 
estimate of a. Where the curved part intersects the value for @ 


Cc 
determines T . 
c A Oe 
The quantity desired is the reduced attenuation @ = a 2 the 
N 


ratio of attenuation due to conduction electrons alone in the super- 
conducting and normal states vs t = T/T» the reduced temperature. 
Subtracting from a. and Op the quantity a. which is attributed to 
impurity scattering, end losses, transducer losses, etc., gives 


a5 =5 Rr 
N T Oo 


Subsequent calculations are very sensitive to the value @ 3? which 
is not easily determined since one must extrapolate asters 40% 
on the reduced temperature scale to t = 0 for zinc. 
The BCS theory gives 
* 
oO = eC TE 
po 
for the variation of ol for the longitudinal attenuation with T, where 


€ (T) is the temperature dependent energy gap between normal and super- 


conducting electrons and k is the Boltzmann constant. What we want is 


asm 





é€ (0). The above expression can be inverted to give 


1 Te 
me °° Cue 


[ 2/0" - 1] 


This is the equation of a straight line passing through the origin 





* 
when 2/a >> 1. Hence, the slope at t = 0 for a plot of y vs t gives 
the reciprocal of the zero degree superconducting energy gap in units 
of kT . 
c 
Data can also be analyzed by comparing the echo height of the 


same echo on each of the pictures. Then for a given n 





A 
a, (db/cm) - eet hn ( ra ) + Oa: 
ni 


One must have some reference attenuation here denoted a, to signify 
the attenuation at the lowest temperature reached. The attenuation can 
be most accurately determined at this point since the maximum number of 


echoes is present. 
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APPENDIX ITI 
PROGRAM RESATT 


The raw data; i.e., pictures, must be first analyzed in the manner 
outlined in Appendix I. Then the computer can be utilized to compute 
the reduced attenuation and the function Y, used to determine the zero 
degree superconducting energy gap, and plot these values separately as 
a function of the reduced temperature. Program RESATT, written in 
Fortran 60 language and compatible with the 1604 CDC computer, will 
accomplish these calculation. The primary purpose of this program is 
to determine the residual attenuation and to find the zero degree energy 
gap. The residual attenuation will be that value which when subtracted 
from the experimental values of attenuation will make Y, for t less than 
0.4, plot as a straight line intercepting the origin. 

For each set of data the program provides for the point plotting 
of computations based on three chosen values of the reduced attenuation, 
each set plotted on a separate graph. On each of these plots a theo- 
retical curve based on the BCS expression, using a chosen zero degree 
energy gap, can be graphed. This is computed from values of €(T)/ € (0) 


vs t, the ratio of the temperature dependent energy gap to the zero 


degree energy gap versus reduced temperature as Rolle 
* 2 
~ = “CCI)EGAP 
Ze 


Yi aera 
Presently 44 values of C(I) are used. More can be computed from the 
reference. 
The title, which must be changed for each set of data, identifies 
the experimenter, date on which datawere taken, frequency, and computer 
Project number. These are within the program itself. The program notes 


follow. 


Code 

TTR(1) - Reduced temperature, t, for theoretical BCS plot 
¥ 

ATS(I) - Reduced attenuation, Q@ , for theoretical BCS plot 
3G 





TY(I) - nec oyaan for theoretical BCS plot 

M - Number of data points for theoretical BCS plot, maximum number, 100 
EGAP - Zero degree energy gap, 2 E (0), (no kT) 

C(I) - E€(TTR(I)/ € (0) 


N - Number of experimental values of T & a, (up to 60) 
TC - Experimental value of Transition Temperature 
AN ~ Experimental value of Q =@Q@ 

n Te 


AZ(O) - Three experimental values of Or 
T(I) - N experimental values of temperature 
ACI) - N experimental values of attenuation 


RT(1) - N values of reduced temperature 


AS (1,1) N values of for AZ(0) - small 
AS (2,1) a _ A - A AZ(0) - medium 
AS(3,1) CG AZ(0) - large 
n O 
AY(1, 1) N values of for AZ(O) - small 
AY(2, 1) 1 AZ(0) - medium 
AY(3,1) y= AZ(O) - large 
fn) =| 
-1 
(O 


Plots: Two sets of point plots 


* 
First: Three plots of Q@ vst, 


Experimental Points AS(1,1I) vs RT(1) for AZ(0) small - CJ 
Experimental Points AS(2,1) vs RT(1) for AZ(0O) medium - © 
Experimental Points AS(3,1) vs RT(I) for AX(O) large - JA 

Each will have the theoretical BCS curve ATS(I1) vs TTR(1) plotted 
on the same graph. 

Second: Three plots of Y vs t, 

Experimental Points AY(1,1) vs RT(I) for AZ(0) small - 
Experimental Points AY(2,1I) vs RT(1) for AZ(0) medium - @ 
Experimental Points AY(3,1) vs RT(I1) for AZ(O) large- A 

Each of these will have the theoretical BCS curve TY(1) vs TTR(I) 
plotted on the same graph. 

mee ee 





For each run a total of 3 + M+ N cards are needed in the following 
order. The last digit of each imput in each column in the stated order 
is to be in the 10, 20, 30, etc., column of the data card. 


Number of cards Imput 
1 M, EGAP 
1 -M Card No., TIR(I), C(I) 
1 N, AN, TC 
1 -N Card No., T(I), ACI) 
1 AZ(1), AZ(2), AZ(3) 


Following is a complete print out of RESATT with the card deck 
composed for the data of 13 Nov 64. Card 000-197 comprise the program 
and cards 200-243 are the Mifhlschlegel values. As outlined, the cards 
to be changed for a different set of data would be 57-61, 199, and 
244-265. 
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APPENDIX III 
PROGRAM VARIGAP 


Once the best value of the residual attenuation and a determination 
of the zero degree energy gap has been made from RESATT, VARIGAP will 
compute and plot the theoretical BCS curves of reduced attenuation vs. 
reduced temperature and the function Y vs reduced temperature for three 
different zero degree energy gaps. On each of these will be the cor- 
responding experimental points for one set of data and the chosen re- 
Sidual attenuation. To the extent that they differ from RESATT, the 


program notes for VARIGAP follow. 


CODE 
EGAP(I) - Three values of the zero degree energy gap (no kTc) 


AS - One value of the residual attenuation 

ATS(1,1) M theoretical values of EGAP - small 
Ado (2,1) reduced attenuation from for EGAP - medium 
ATS (3,1) the BCS expression EGAP - large 
yc, L) M theoretical values of EGAP - small 
meee, L) Y from the BCS expression for EGAP - medium 
ec, L) EGAP - large 
AS (J) - N experimental values of reduced attenuation 

AY (J) - N experimental values of Y 


PLOTS: Two sets of curve plots 
First: Theoretical reduced attenuation vs reduced temperature 


ATS(1,1) vs TTR(I) for EGAP - small - LJ 
ATS(2,1) vs TTR(1) for EGAP - medium - © 
ATS(3,1I) vs TTR(I) for EGAP - large - A 
Each of these will have point plotted on them the N experimental 


values of AS(J) vs T(J). 


Second: Theoretical Function Y vs Reduced Temperature 


TY(1,1) vs TTR(L) for EGAP small 
TY(2,1) vs TIR(1) for EGAP medium 
TY(3,1) vs TTR(I) for EGAP large 


Ai5- 





Each of these will have point plotted on them the N experimental 


values AY(J) vs T(J). 


For each run a total of 3 + M + N cards are needed plus the correct 


title cards. The last digit of each imput in each column in the stated 


order is to be in the 10, 20, 30, etc., column of the data card. 


Number of Cards 


Imput 
M, EGAP(low), EGAP(medium) , EGAP(large) 
Gard, No., TIRCL), CG) 


N, AN, TC 
Card No., T(J), ACJ) 
AZ 


Following is a complete print out of VARIGAP with the card deck 


composed for the data of 13 Nov 64. Cards 000 - 188 comprise the program. 


As outlined, the cards to be changed for a different set of data would 


be 57 - 61, 190, and 235 - 257. 


See 





09€000 
0Sse000 
07€ 000 
o0¢e000 
0z¢€000 
OTEO00 
00¢€ 000 
062000 
082000 
012000 
092000 
0S72000 
072000 
O0¢z000 
022000 
OTZ000 
002000 
06TO000 
O8TOOO 
OLTOOO 
09TO000 
OSTOOO 
07TO0O0 
O€TOOO 
O02TO00 
OTTOOO 
O0OTOOO 
060000 
080000 
020000 
090000 
0S0000 
070000 
0£€0000 
020000 
OTOOO00 
000000 


ZS ISS IS((1 ULL SI 
WST=I IT Od 
€ l= Pleo 

(O°OT4S) LVWYOS 
Zv*OT avau 
(0°OTSZ*OTI) LYWYOS 
(CTIVECTILS1)%6 GAY 
N‘T=r 8 0G 
DLONVONSS GVAY 
(CLIDSCTIMLLEII 66 GVaY 
W*T=C L Od 
(O°OTSESOTI) LVWYOS 
dVOS*WS0S avay 
°0=(1)dv93 

€*T=I 9 Od 

SO=(P)AV 

"O=(F)SV 

*Os(7)1¥ 

O€*I=F S OG 

*O=(1)V 

‘o=(1)1 

O9*I=1 4% Od 
"O=(ISPIAL 
*O=(I SF )S1V 

e*t=f € Od 

*O=(1)D 

eO=(1)YLI 

OOT*T=I € Od 
H@=(1)SLILI 

21*I=I Zz 0a 

SO=( TDA 

*O=(1)X 

006*I=I T Od 


(EVdVOSS(OEILYS(OLCIAVS(OEISVS(O9DIVS (O09) LS(OOTST 
CMALS(OOTSEISIVIE(OOTIDS (COOTIYLL¢S (CTISILILI£(006)AS(006)X NOISNIWIQ 


dVOTYVA WVeD0dd 


OS 


ZINOO* JSTeec odor ®* — 


eet 


By AG | 





0€2000 
021000 
OTL000 
002000 
069000 
089000 
019000 
099000 
0S9000 
079000 


0€9000 


029000 
019000 
009000 
06000 
085000 
015000 
09¢000 
0S¢S000 


0475000 


0€S000 
02S000 
OTSO00 
00S000 
0647000 
087000 
017000 
0947000 
0¢%7000 
0477000 
0€%7000 
027000 
0147000 
007000 
06¢£000 
08¢000 
OL€000 


O=XVX3IGOW 
O=LHOTYAI 


O=dNXI. 


Z°O=31VDSA 
Z*O0=3 1VDSX3 
JYUNLVHS=(Z2TIZILILI 


YAdW3SL GHB=(TTISAILILI 
JONGSY SHB=(OTISZILILI 


A NOTLVNH8=(6)39TLIL1 
N3LLV G3SHB=(8)3STLILI 
ONGSAY H8=(LZ)37L111 
SOHB=(9ISATILILI 

WOE 79 HB=(SIZSILILI 
AON €T TH8=(¥7)3SITLILI 
c£O LOSFHB=(CISILILI 
Odd ZINHB=(c)SILILI 
09 HB=(TISZILILI 
H8=139V1 


SANILNOD 
CANN) S9OOT/* T=( PD AV 


91 OL 09° 


°O=(F)AYV 
GI*7t*7l(*T-ANNG) SI 
*T-(FISV/°?C=ANNG 


(ZV-NVIS(ZV-(fF IVI =RHCC OSV 


Dis(f Ler yiy 
N*T=f 9T Od 

Oc=N 
CLSeTtT Se T(OE-N) dl 
JINNILNOD 


(INNG) SOOT PE THC ISP dAL 


7S *€S%ES (°T-ANNG) JI 
*T-(CISFISLV/*°2=ANNG 


C*TH( CCD YLLle®* 2S (TL IDe( fF) dV9OIIAdX3S)/° CHIT SP SLV 


ttl ol O9 
°O=( ISP IAL 
*O=(1 Sf )S1LV 


oT 
ST 


yT 


Et 
cl 


Tt 
7S 
cS 


Ss 
TS 





OTTTOO 
OOTTOO 
O060TOO 
O8s0TOO 
OLOTOO 
090TOO 
OSOTOO 
0vVOTOO 
OcOTOO 
OZOTOO 
OTOTOO 


000 TOO 


066000 
086000 
016000 
096000 
0S$6000 
076000 
0¢€6000 
026000 
O0T6000 
006000 
068000 
088000 
0L8000 
098000 
0S8000 
078000 
0€8000 
028000 
OT8000 
008000 
062000 
082000 
OLL000 
091000 
0SL000 
0471000 


W°T=f T2 Od 
T-W=S1dWNN 
€*T=I 42 0G 
CT=HOITHI 


JYNLVYHB=(ZTISILILI 
JdW3Ll GSH8=(TTISILILI 
INGSY SAHB=(OTISILILI 

((T-sxVHH8=(6)3 LI 11 

d1V/c )OOHB=(B8ISILILI 
VW/THA HO=CLISAILILI 


°O=(F IA O2 


*O=(f)X 

N*Tt=f O2 Od 
(LSVISGTYOT SHOTHI SSGIMI * XVASGOWS XVXAGOWS LHOTYAT SdAXxIT 
SITVISASSAIVISXASAILILI GS VAGV 1S AdALISAYNDIGOWS ASXSSILGWANIMVYG WIV 
C+ =3dAlLI 

€=AYNDAGOW 

(FISW=(FIA 

(CT) LY=(F)X 

N*T=F 6T Od 

N=S1ldWNN 

°O=(T)A 

*"O=(T)X 

W*T=f 8T Od 
(LSVISQT YO! SHOTHI SSGIMI § XVASGOW * XVXSGOWS LHOTYAT SdAXIT 
SJIVISASAIIVISXSSSILILIS VSA9V 1S AdALISAYUNDGOWSASXSSLdGWANIMVYG 11V)D 
O=3dAlI 

T=AYNDIGOW 

(FST ISLV=EICIA 

(TF YLL=a(C KX 

W*T=fF LT Od 

W=S1dWNN 

e*t=I O¢ Od 

O=1SV1 

T=qlyo] 

8=HOTHI 

8=30IMI 

0O=XVAAIGOW 


6T . 


ST 


LT 


mmr meme cm ts rn ES en 


~49- 





OL%7TOO 
09%7T00 
0S7T00 
Ov4TO0O 
O€+TO0 
O027T00 
OT7TOO 
0047TO0 
O6€ TOO 
Ose tod 
OLE TOO 
O9€T00 
OSseTOO 
OvE TOO 
OeceTOO 
OceTOO 
OTETOO 
OOE€TOO 
062TOO 
O8cIo00 
OLZTOO 
O097TOO 
OScT00 
O72 TOO 
O€cTOO0 
O2zTo00 
OTZTOO 
002TOO 
O6TTOO 
O8TTIOO 
OLTTOO 
O9TTOO 
OSTTOO 
OFT TOO 
O€TTIOO 
OcTTOO 


\ 


(NS THIS (CTOSVECTILYSCI)VS CTI L)¢ce INIUd 
(//%WHd1V JSYUNLVYSsdwW3sl VT 
Hd 1V SUNLVYSdW3L | H98) LVWYOd TE 
Te LNIdd 
(ag35Nd034y i 
HZl) LVWYOS OF 
O¢€ LNIUd 
(//S°84*= ZV ONV HITSS°8se= T 
SL * H6*S*°8s%= NY HO) LVWYHOS 62 
ZVSDLSNVS6C INI Ud 
(YO4 SSANIVA WINSAT 
WIY¥3SdX3 | HITZ) LVWYOS 82 
: 8c INIUd 
qGQaadnat 
TH88) LVWYOsSO0Le 
Le iNI&d 
"O=(CIA V2 
"O=()X 
N*t=f %¢ Od 
(LSVISGQIYST SHOFHI § S0IMI * XVASZGOW S XVX3IGOWS LHOTYAT Sd AXIT 
SSIVISAS A IWOSXASAIILILI SS VAGV IS SdALISAYUNDGOWSASXSSLGWANIMVYDO T1V> 
€+T=3dAlLI 
€=AYNIIGOW 
(TIAVEITIA €2 
(T)LY=(F 9X 
N*T=f ¢€2 Od 
N=S1LdWNN 
"O=(F IA C2 
°O=(fF)X 
WST=f ¢2 Od 
(LISVISACIYOT SHOTHI SS0IMI Ss XVASGOWS XVXSIGOWS LHOTYAT SdNXIT 
SJTIVISASSIIVISXSSSIILILIS T99V1S SdALISAYNDIOWSASXSSLGWANIMVYD W1V5 
O=3dAlLlI 
T=AYNIGOW 
(CSTIALHAIC IA 2 
(rc yali=a (er dx 


(J///9N YNLVASBSdgWSBl 
34 S A * VHd 1/9 


-50- 





Ov8T00 dV93*8E LNIUd 


O€8TOO LE ANIUd 
OZ8T00 9€ LNIUd 
OT8TOO : (AL I 
0O08TOO JYUNLVYSdW3aL HOL) LVYWYOS 2» 
O6LT00 ( 247 INIUYd 
O8LT0O 7€ INIUd 
OLLTOO €€ iINIYd 
O9LTOO (N° THIS (IDAVS CTI LYS CT VEC II L)°ZE LNIUd 
OSLTOO (//4V JYNLVYIdWaL vt 
O¥LT00 Hdtv ~2YNLVYSdWaL Hg) 1LVWYOS TY 
O€LTOO Tv LINIUd 
OZLTOO O€ INIUd 
OTZLTOO , ZVSDLSNVS62 LNIUd 
OOLTOO | 8z LNIUd 
069T00 (JJV/IUYNLVAIAdGwWIAL q3T 
089100 PNAAZRBY S A A ; ITHE8) LVYWHOS OF 
0L9TOO Ov INIUd 
099T00 (9°E€TAHZS9*LTH49°E74) LYWYOS 6E 
0S9TOO (WS T=ISCISENSLYVS CISZ)SLVECISTISLVSC I) ULL) S6€ LNIUd 
0”9T0OO (//€°ETAZSE*°6S4) LVWHOS BE 
O€9TOO dv93*8e€ ILNIYd 
0z9T00 ye "ee _=dvot 
OI9TO0O 3 HG8) LYWYHOS LE 
009TOO Le ANIUYd 
06ST00 ( satasheinsinateataeiamneatal atasieeieataaiaaieniennnant T 
O8st0oOo -——-- HL8) LVWYOS 9E 
0OLSTOO 9€ INIUd 
09STOO (xVHd1V I 
0SSTOO JUNLVYSdWIAL HZL) LVWYOS SE 
OrSTOO gee Se INIUd 
O€STO0-- =~ ~  (@3aDNna3a4y Hev) LWWYOS vE 
~ 026100 ve INIUYd 
OTSTOO (S//SANIVA WOIL3T 
0OOSTOO %sO3H1L H89) LVWYOS EE 
06%T0OO €€ LNIYd 


087100 (///(9°7T3SS°9TASS*°CTASS°244)) LVWHOS ct 





002zZ00 6808 °0 ZL°O O€ 
061Z00 88Z78°0 0L°0 62 
O8TZ00 . TL+9°0 g9°0 82 
OLTZ00 07998 °0 99°0 L2 
09TZ00 — ss Soe _9628°0  79°0O 9 
0Stz00 6€68°0 29°O S2 
OT zZ00 0L06°0 09°0 42 
0€1zZ00 0616°0 8S5°0 €z 
OZ1Z00 6626°0 9$°0 22 
OTTZ00 66€6°0 4S *0 tz 
OOLZ00. -. 88476 °O Z25°0 02 
060200 69S6°0 0$°0 6T 
080Z00 . 1796°0 84 °O eI 
010200 70L6°0 94 °0 LT 
090200 0926°0 77° 91 
0S0200 6086°0 Z4°O ar 
040200 0S86°0 O%°0 oT 
0€0Z00 6886 °0 8€°0 El 
020Z00 G166°0 9€°0 ral 
010200 , 8E66°0 vE°0 IT 
000zZ00 1$66°0 Z€°O OT 
066100 1266 °0 0€°0 6 
086100 Z866°0 gZ2°0 8 
0L61T00 6866 °0 92°0 L 
096100 7666 °0 ¥Z°0 9 
0S6100 1666°0 ZZ2°O S 
046100 6666 °0 02°0 4 
O0€6T00 0000°T 81°0 € 
026100 0000°T 9T°0 Z 
016100 0000°T 71 °0 I 
006100 89°E QS°¢ Sy °E be 
068100 

OsgI00 Gna 
OL8T0O aNa 
098100 dOLs 


O0S8T00 (WSTHIS CI SEVALSCISEZSIALS CIS TIALS (I) Ylid*sé6e INIud 


2 we re eee eee ee 


Boe 





ee 


096200 
069200 
047G¢ 200 
0€S200 
029200 
O0TS200 
00S200 
067200 
0847200 
OL7200 
094200 
0S%7200 
0747200 
0£%7c00 
027200 
OT7200 
007200 
06¢200 
08¢€200 
OL€200 
09€200 
0S€200 
O0v7€Z200 
O€€200 
O0ze200 
OTEzo0 
00¢€z200 
062200 
082200 
042200 
092200 
062200 
072200 
0€72c00 
O2cc00 


OT2c00 


c60°%E 
c60°E€ 
cOO0°%E 
C60°E 
S90°E 
768°C 
SZG we 
O8c°%e 
Sv76°T 
799°T 
Lev*t 
e8e°o 
G77°O 
777 °O 
9¢7°0 
7e7°O 
7lL7°O 
T67°0O 
7747 °O 
7l1S°0O 
LT8°O 
0000°0 
9¢72°O 
9TvE°O 
87T7°O 
671%7°O 
€92S°0 
STLS°O 
LTS 
08479°0 
OT89°0O 
OTTL°O 
98€L°O 
0791°O 
718L°O 


TOT 
61T6°0 
898°0 
€£8°O 
eTs°*o 
c08°0 
€Ll°O 
€91°0O 
7EL?O 
€0l°O 
€89°0O 
O€€ °O 
6€€ °O 
Ele °O 
8LE°O 
T8¢e°O 
98¢€°0O 
O6€ °0 
96€ °O 
cc7°O 
c60°E 
00°T 
86°0O 
96°0 
76°0 
c6°O 
06°0 
88°O 
98°0 
78°0O 
c8°O 
08°0O 
8L°O 
91°O 
7L°0O 


IANO TMH OM O©OAO 
ae | 


Oct°O 


02 
6T 
ST 
LT 
9T 
ST 
wT 
eT 
cl 
TT 


O¢ 


ty 
CY 
TY 
he 
6€ 
BE 
Le 
9¢€ 
Se 
7€ 
ce 
ce 
Te 


i> 


aia 





APPENDIX IV 
PROGRAM ZEROGAP 
This program provides for a least square fit of experimental points 


to the BCS function 


Ss 2 
a €(T) 
N 
kT 
AL he | 


by minimizing the square of the vertical distances from the curve to 





the experimental points. This is done using two parameters; the re- 
Sidual attenuation and the zero degree energy gap. The output of the 
program is the value of these two parameters. Since in the BCS ex- 


pression Or. and @. involve only the attenuation due to electron-phonon 


N 
interaction, all other causes of attenuation must be subtracted out. 


Thus, 
| 
0 | ae 
AO, ee 
L + 1 
or 9 
1 22 
a ea (lee) €(T) 
L kT + | 


where Q' and On would be on a BCS curve shifted an amount Ot due to 
residual attenuation. 
What is desired is to minimize the residue 
N 
Re » (a-a')* 
in 
where Q@ is an experimental measurement of attenuation and N is the number 
of experimental points. 
Let the following quantities be defined which will be used in the 
program. 
Q - ACI) - Experimental measurement of attenuation 


t - RT(1I)- The reduced temperature corresponding to the temperature T(1I) 
at which A(I) was measured. 


Eis 





a AZ - Residual attenuation - one of two parameters to be determined 


C(I)-US(I) - The Muhlschlegel values of €(T)/ €(0) vs t. 


By use of an interpolating program, the values of C(I) were obtained 
from the Mtihlischlegel values for 0.191 <¢ t € 1.000 in increments of 
0.001 and incorporated in the program. In the printout of the program 
following they are cards 200 - 910, where for example, 200 means 
t=0.200 and the number in the 10th column is C(I) for the t. 
Gq) /t = UC(I) 


The equation for @' can then be written 


a' =a Qa -a ) Se 
mato ( No C(I): 2 €(0) 
pre + 1 
Now let FER(I) = fe) ee, ee eee 
C(I)- 2 €(0) UC(I) «EGAP 
at EXP | ——4——=— | + 1 
L + 1 2 


where EGAP is the zero degree energy gap = 2 €(0) - the other parameter 
to be determined. Then the residue becomes 
) 
Kee y [ Ac1)-az-(avrazy2PeR(1 | : 
{=1 


From this an explicit expression for AZ can be gotten from the 


condition 
aR 


DAZ = QO. 


Carrying out the steps gives 


N 
21 | (2ANFER(I)-A(T) | | 2FER(I)-1| 
AZ = ~ 


=, [ arercry-1] [ 2rerqy-1] 


in which terms cannot be canceled since they are under the summation 


signs. In this expression only FER depends on EGAP. 
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AR 


The other condition for minimizing the residue is @ EGAR 


Carrying out the steps gives 


OR 
a EGAP 2 2. 


af 


N 
| A(1) -AZ-2(AN-AZ)FER(I) 


x }2(AZ-AN) = FER(FER-1) | 


So what is needed is to minimize Y. [ (2ANPER(T) -A(I) ) -AZ(2FER-1) | UC 
x | FER(FER-1) | [ AN-AZ | 


subject to the above value of AZ. The remainder of the program notes 


are the following; 


ONE (1) 
TWO(I) 


ZFER(I) - 1 
2ANFER(I) - ACI) 


The imputs needed are the best guess of EGAP; AN, the experimental 
determination of the attenuation in the normal state; TC, the transition 
temperature; and the N data points of A(I) vs TCI). 

One card contains N, AN, TC, and EGAP, the last digit of each is in 
column 10, 20, 30 and 40 respectively. Next follow the N data cards of 
ACI) vs TCI); one card for each set with the last digit of the card number 
in column 10, that of T(I) in 20, and that of A(I) in 30. A print out of 
the program follows at the end of the appendix. The imput cards follow 
the blank card which must be added to the end of the program cards. 

Table III gives the results of program Zerogap and compares them 
to those of program Resatt. Figs. 15-26 give the plots of the BCS curve, 
using the values of residual attenuation and the energy gap determined 
from Zerogap and the point plots of the corresponding experimental data. 

It is evident that the program does what it is designed to do. 
However, the values of residual attenuation are markedly different in the 
30 mcs/sec data from that obtained by extrapolation of the curve of 
attenuation vs temperature. For the 30 mcs/sec, 5 Jan 65, semi-log data 
the value obtained is unrealistic in that it is greater than most of the 
data points. Only the attainment of lower values of reduced temperature, 


=S6- 





t, can resolve the question as to whether or not the value of residual 
attenuation obtained from Zerogap is realistic. 

The program seems to work best when there are data points over 
the whole range of t; e.g., 30 Nov 64 and 5 Jan 65, 10 mcs/sec data. 
Also, when the points are not too scattered as in the 30 Nov 64 data, 
the results are reasonable. The values of (2.85 and 2.82)kT obtained 
for the energy gap for these data are reasonable bearing in mind the 
previous observation that for higher t the experimental points lie above 
the 3.50 kT. BCS curve. As one can see the lower the value of the energy 
gap, the higher the BCS curve. 

This program can be of use in interpreting data if one realizes 
that the analysis depends completely on the BCS theory and the simplifying 
assumptions used in obtaining the expression for the reduced attenuation 


and masks the possibility of other phenomenon taking place. 
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is 
13 


Ww 


Nov 
Nov 


Jan 
Jan 


Jan 
Jan 


64 
64 


65 
65 


65 
65 


TABLE III 


COMPARISON OF ZEROGAP AND RESATT RESULTS 


FREQUENCY 
Mcs/sec 


30 
30 


30 
30 


10 
10 


TYPE 


Semi-log 
Echo Heights 


Semi-log 
Echo Heights 


Semi- log 
Echo Heights 


ZEROGAP 
RES IDUAL ENERGY 
ATTENUATION GAP 
db/cm 2€(0)/kT. 
Ori 2.8544 
0.1507 2 ones | 
0.8755 7.2968 
0.4387 PAS SOY 
0.7166 3.5941 
0.6097 Zelvst 
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RESATT 


RES IDUAL 


ATTENUATION 


db/cm 


Oo © 


-320 
. 300 


D005 
25985 


. /00 
lS 


ENERGY 
GAP 
2€(0) /kT. 





- Te, 


Semi-log Blot Data , 
a, = 0.1772 db/em a 
2€(0) = 2./8544 AT. 
& 
© 
i i 
cy 
a} 
P. 


619 


ns. 


: 
‘ 
: 


N“OCALE = 2,08E-@1 UNITS”INCH. 
"SCALE = 2005-01 UNITSVINCH, 


GONCZ . PROJECT 8321 13 NOU 64 3QMCS 
REDUCED ATTENUATION YS REDUCED TEMPERATURE 


Fig. 15 
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Semi-log Plot Data | 
Ot = 0.1772 db/cm 
2€(0) = 2.6544 kT. 


meeceLe = ZOOF- O1 LINTT S-“INCH. 
Y-SCALE = 2,00E-01 UNITS INCH. 


GONCZ PROJECT 8321 18 NOUV 64 SBMCS 
f= IA(LOG(Z’“ALPHA*- 1) YS REDUCED TEMPERATURE 
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Fig. 16 
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P“SCPLE = ZOOE-O1 UNITS INCH. 
-SCALE = 2,00F-01 UNITS” INCH, 


GONCZ PROJECT 8321 153 NOU 64 3QMCS. 
REDUCED ATTENUATION YS REDUCED TEMPERATURE 


| Fig. 17 
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Echo Height Data 
oe = 0.1501 db/cm 
2€(0) = 2.B17 KT. 


NSCALE = 2,.00E-@1 UNTTS”INCH. 
Y-SCALE = 2.Q0F-01 UNITS“INCH, 


GONCZ PROJECT 8321 13 NOU 64 3EMCS 
f= 1A(LOG(2“ALPHA*- 1) YS REDUCED TEMPERATURE 
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Fig. 18 
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Semi-log Plot Data 
cae = 0.8755 db/em 
2€(0) = 7.297 RT, 
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K-SCALE = 2.005-91 UNTTOVINCH, 
Y-SCALE = 2.00F-@1 LUINITS/INCH. 


GONCZ PROJECT @321 5 JAN 65 3QMCS 
REDUCED ATTENUATION US REDUCED TEMPERATLE: 


Fig. 19 
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Le) : 
Semi-log Plot Data 
a. = 0.8755 db/cem 
ZECO) = 7/2297 RT. 
ra) 
©! 
a) 
| feo (| | 
we 086 =~ 828 910 | 


A-SCALE = 2.082-01 UNITS-“INCH. 
feoereE = 2.082£-01 UNITS“INCH. 


GONCZ PROJECT 98321 5 JAN 65 32MCS 
f= VLOG /ALPHA*- 1) YS REDUCED TEMPERATURE 


Fig. 20 
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Echo Height Data 

a. = 0.4387 db/cm 

2€(0) = 2.530 KT. 
S : 
rT | 7 
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G39 
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CO2 
A-SCALE = 2,00F~-@1 UNITS- INCH, 
_ f-SCALE = 2,002-91 UNITS. INCH, 


GONCZ PROJECT BS21 5 JAN 65 32MCS 
REDUCED ATTENUATION YS REDUCED TEMPERATURE 


meee 21 
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Echo Height Data 
oe = 0.4387 db/cm 
2€(0) = 2.530 RT. 
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,OCRLE = 2008-81 UNITS. INCH, 
Y°SCALE = 2,00F-01 UNITS. INCH, 


GONCZ PROJECT OS21 5 JAN SD 3EMCS 
Y=1-(LOG(2’ALPHA*- 1) YS REDUCED TEMPERATURE 
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Semi-log Pilot Data 
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“SCALE = Z.O0E-O1 LINITS. INCH. 
SCALE = 2.08E-O01 UNITS INCH, 


GONCZ PROJECT 8321 5 JAN 65 1QMCS 
REDUCED ATTENUATION YS REDUCED TEMPERATURE 
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Semi-log Plot Data 
oe" =O leO db/cm 
2€(0) = 3.$94 kr. 
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NOOLALE = 2,00F-01 UNITSINCH. 
Y-SCALE = 2022-31 UNITSVINCH. 


GONCZ PROJECT 0321 5 JAN 65 10MCS 
f= 14 LOGW2“ALPHA*- 1) YS REDUCED TEMPERATURE 


Fig. 24 
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Echo Height Data 
Ot = 0.60917 db/cm 
2 €(0) = 2.1118 kT. 
ran 
© 
Pocomee = %O0E-01 UNTTS-INCH 


6 A 
Sees = 2.00--01 UNITS INCH. 


GONCZ PROJECT 8321 5 JAN 65 1QMCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 





Figs 725 
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—  GONCZ PROJECT @321 5 JAN 65 1@MCS 
Y=1-(LOG(Q/ALPHS*-1) US REDUCED TEMPERATURE 
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APPENDIX V 
PROGRAM AVERGAP 


This program provides for a least squares fit to the BCS function 
using one parameter; the zero degree energy gap with the residual atten- 
uation specified. Thus the program and the program notes are the same as 
those of Program Zero-gap, except that the residual attenuation is part 
of the imput data rather than being derived from an analytic expression. 
The program will compute the zero degree energy gap which minimizes the 
square of the vertical distance from the experimental points to the BCS 
curve for three different values of residual attenuation. Thus, the 
expression 


, | (2ANFER (I) - A(I)) - AZ(2FER(I) -1)| UC | FER(FER-1]| AN-Az| 
‘=I 


is minimized for each given AZ. 

The type of data cards, their format, and order of imput are the same 
as for Program Zerogap, with the exception that one card need be added. 
This is the last card and contains the three values of residual atten- 
uation; the last digit of each to end in column 10, 20, and 30 respectively. 
These values are to be determined as in Program Resatt; i.e., from an ex- 
trapolation to T = 0 on a plot of the experimental data of attenuation vs 
temperature. 

Table IV gives the results of Program Avergap and also a summary of 
all the computer calculations. It is seen that the results of Avergap 
are more consistent with those of Resatt than are those of Zerogap and 
again reflect the fact that the experimental points for higher t lie above 
the BCS curves used in Resatt. The results for the data of 5 Jan 65, 30 
Mcs/sec still seem unreliable with such a high energy gap. This is due to 


the fact that there are data points for low t only. 
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